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Abstract 

Using holographic renormalization, we study correlation functions throughout a 
renormalization group flow between two-dimensional superconformal field theories. 
The ultraviolet theory is an N = (4, 4) CFT which can be thought of as a sym- 
metric product of U{2) super WZW models. It is perturbed by a relevant operator 
which preserves one-quarter supersymmetry and drives the theory to an infrared 
fixed point. We compute correlators of the stress-energy tensor and of the relevant 
operators dual to supergravity scalars. Using the former, we put together Zamolod- 
chikov's C function, and contrast it with proposals for a holographic C function. In 
passing, we address and resolve two puzzles also found in the case of five-dimensional 
bulk supergravity. 
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1 Introduction 



In the AdS/CFT correspondence, renormalization group (RG) flows of a <i-dimensional 
conformal field theory are described by domain wall solutions of the dual (d + ^-dimen- 
sional bulk supergravity theory, see e.g. PQ and references therein. Physically, the domain 
wall solution can be thought of as a shell of matter, where the metric becomes asymptot- 
ically AdS far from the shell. In terms of limiting behavior of the dual field theory, the 
asymptotic AdS length scale gives the central charge at the ultraviolet conformal fixed 
point of the field theory. Should the field theory be conformal also in the infrared, the 
bulk space is asymptotically AdS also in the deep interior of the shell, with the inner 
AdS length related to the infrared central charge of the field theory. In the supergravity 
scalar target space this means that the flow does not run off to infinity but goes down 
to a minimum of the potential and stops. The solution is nonsingular everywhere, so for 
weak fluctuations and small curvature, the supergravity approximation may be trusted 
throughout the flow. 

The previously studied five-dimensional examples of such conformal-to-conformal flows 
could only be given numerically 0. This has remained an unsurmountable obstacle for 
further application of holographic renormalization methods, in particular for the compu- 
tation of correlators in nonsingular flows. The few flows that are known exactly, on the 
other hand, all run off to infinity in the scalar target space, which creates singularities 
at finite distance in the bulk spacetime Ej. On the field theory side, these solutions 
describe theories that confine in the infrared. In [3] , for three-dimensional bulk, we found 
the presently only known exact solution describing an RG flow to an infrared fixed point. 
The present paper is about correlators in this smooth RG flow between conformal fixed 
points. Our main result is the computation of two-point functions of the stress-energy 
tensor and operators dual to the supergravity scalars. 1 

To be specific, the flow we study is a solution of the three-dimensional SO (4) x SO (A) 
gauged supergravity |H| with equal coupling constants and 16 supercharges. The dual field 
theory is a large (or "double") N = (4,4) superconformal theory with SU(2) A current 
algebra and equal levels of the two SU(2) 2 factors. It has an alternative realization as 
a symmetric product of U(2), N — 1 super WZW models [7], and in terms of branes, 
it is the worldvolume theory on the intersection (along the Dl-branes) of two D1-D5 
systems with equal D5-brane charge. This brane setup has near-horizon limit Ad S3 x 
S 3 x S 3 x S 1 , where one can recognize the two SO (4)'s as transverse rotations of the 
two brane systems. Finally, the RG flow we describe is driven by a relevant operator 
O q of ultraviolet dimension 3/2, that breaks conformal symmetry and leaves 1/4 of the 
supersymmetry. 

1 Guidelines for selected reading of this paper are given at the end of the introduction. 
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Let us briefly expand on our motivations for this work. On a purely methodological 
level, we want to explore the extent to which three-dimensional gauged supergravity can be 
used to describe deformations of interesting but complicated two-dimensional CFTs, such 
as worldvolume theories on intersecting D-branes. In the end, we hope to have provided 
an example that three-dimensional gauged supergravity does provide a powerful addition 
to the arsenal of methods in two-dimensional CFT. In particular, on the supergravity 
side we are able to compute quantities that are very difficult to compute in the deformed 
CFT, e.g. (TT) correlators and the C function. It is intriguing that on the field theory 
side, those computations seem to be difficult but not impossible, which makes our results 
predictions for supersymmetric (and perhaps integrable) deformations of two-dimensional 
CFT. In addition to this motivation, our flow may be viewed as a toy model in which 
to address some questions raised in flows of higher dimensional theories, such as four- 
dimensional N = 1 super Yang-Mills theory. The toy model idea proves to be useful, 
since we do find some new insight into old mysteries encountered in holography with 
five-dimensional bulk supergravity (see section EJ. 

For completeness, let us mention that in [S] we did not just find the supersymmetric 
flow but we also found a stable nonsupersymmetric fixed point of this theory (in fact, 
the first stable nonsupersymmetric fixed point; later on several others were found in the 
three-dimensional maximal theory jHJEl)- A flow to such a nonsupersymmetric fixed point 
may be relevant to some aspects of black-hole physics, and this was part of the motivation 
of 5J. However, although we hope to discuss black- hole physics in future work, here we 
focus entirely on the supersymmetric flow and the complementary motivations mentioned 
in the previous paragraph. 

To obtain the correlators, we use the formalism of holographic renormalization fUJ 
El El El- It is a framework to reliably compute correlators along RG flows of quantum 
field theories using the AdS/CFT correspondence. In particular, it allows one to compute 
one-point functions in the presence of sources, to check that they obey the requisite 
Ward identities throughout the flow, and to compute power-law terms which were usually 
dropped in the old prescription E] • I n a theory that is conformal also in the infrared, 
such as the one considered here, surely a basic requirement of the formalism one wishes 
to apply is that it is sufficiently restrictive to single out the correct asymptotic power- 
law behavior automatically. This is true for holographic renormalization, as we show for 
example in section lo~Il As a consequence of being able to compute one-point functions, one 
can also make sure one is using a renormalization scheme that preserves supersymmetry, 
since (Ty) = is then expected to hold in the background. Finite counterterms have to 
be added to ensure this, corresponding to a selection of a renormalization scheme that 
preserves supersymmetry in the dual field theory. 

After this brief introduction to the formalism, we proceed to summarize the new 
results in this paper. We extend the analysis of [TT1 Ej to also include supermultiplets 
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consisting of "inert" scalars $ (meaning they have vanishing background — as opposed 
to the "active" scalar Q which carries the domain wall background). As it turns out, the 
treatment of the inert scalars raises some new conceptual questions, such as the proper 
choice of finite counterterms and the distinction between inert scalars with conformal 
dimensions A + and A_, corresponding to the same mass in supergravity, as we shall 
shortly discuss. Through the coupling to the active scalar Q, the inert scalars change 
mass along the flow to the fixed point. This coupling also allows for (in fact, requires) a 
new type of finite counterterm, of the form Q 2 <& 2 . We compute the coefficients of these 
terms for all inert scalars, as well as the coefficient of the finite counterterm Q 4 (first 
displayed in [11). The terms of type Q 2 $ 2 can be seen as natural generalizations of the 
finite Q A term, but the coefficients of the Q 2 $ 2 terms cannot be determined by evaluation 
on the background, since $ itself vanishes on the background. In section 15.31 we will see 
how to compute those coefficients, using a supersymmetry Ward identity for the two- 
point functions of superpartner inert scalars; the result is listed in (I5.27j) . By adding all 
counterterms and taking the cutoff e back to zero, we compute the renormalized action. 
As usual, all correlators are computed by functional differentiation of this renormalized 
action with respect to the boundary sources, and then setting the sources to zero. 

The next result is more conceptual. It was noticed in j3] that the effective potential 
V appearing in the fluctuation equations of the inert scalars can be expressed in terms 
of a simple prepotential in the sense of supersymmetric quantum mechanics (susy-QM). 2 
This property is very useful, since it guarantees the absence of tachyonic fluctuations. 
Similarly, such prepotentials were found in the known five- dimensional examples [Tj)]: 
the general existence of susy-QM prepotentials seemed in need of further explanation. 
In section 15.11 we show that the existence of susy-QM prepotentials follows from the 
preserved N = 1 supersymmetry of the background flow, and that the prepotentials may 
be directly extracted from the fermionic mass term of the underlying gauged supergravity. 
This argument extends readily to higher dimensions. This is one aspect in which we see 
that the toy model does provide useful information in higher dimensions. 

We then proceed to determine the supersymmetry Ward identity relating two-point 
functions of a pair of superpartner scalars. This is useful as it provides a way to distinguish 
between different conformal dimensions associated to the same supergravity mass. It is 
well known that scalars of mass in a certain range (here — 1 < m 2 < 0) can correspond 
to two different solutions A± for the conformal dimension of the dual operator. For the 
active scalar, the choice of A + or A_ is the difference between whether the background 
describes an operator flow or a vev flow [17J, but this is not so for the inert scalars. 
Our scalars are precisely in this range, and unlike in previously studied cases, there are 

2 In this paper, the notion of prepotential always refers to this supersymmetric quantum mechanics 
function U encoding the potential of a fluctuation equation as V = U'+U 2 , as opposed to the superpotential 
W which describes the background potential V of the active scalar as V — \(W') 2 — 2W 2 . 
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now two representation sectors with the same quantum numbers, so group theory is not 
sufficient to make the distinction. Using the fact that the correlator asymptotics depends 
directly on A and not just on the supergravity mass, and following the correlators from 
one fixed point to the other, we are able to distinguish between A + and A_. 

We then derive the fluctuation equations for the inert scalars, active scalar and met- 
ric around the domain wall solution. In all previously studied cases, those fluctuation 
equations were hypergeometric. Instead, we find that around our domain wall solution all 
fluctuation equations, for inert and active scalars as well as for the metric and vector fields, 
reduce to a slightly more complicated equation, the biconfluent Heun equation, which de- 
scends by confluence from a Fuchsian equation with four regular singularities. We devise 
some methods to solve this equation; the mathematics is relegated to Appendix iBl where 
in particular we point out a simple and efficient way to compute the sought-after coeffi- 
cient numerically. This allows us to achieve our main goal: the computation of two-point 
correlation functions throughout the renormalization group flow. 

The final part of the paper concerns the computation of a C function, i.e. a function 
that is monotonic as a function of RG scale along the flow and interpolates between the 
central charges at the conformal fixed points. The general existence of such a function 
can be very useful to map out the space of field theories, for instance to find well-defined 
universality classes. In two dimensions, Zamolodchikov has given a general construction in 
terms of the two-point correlators of the stress-energy tensor |18j . Since Zamolodchikov's 
proof of monotonicity relies heavily on the lack of distinct tensor structures for stress- 
energy 2-point functions in two dimensions, it has no straightforward generalization to 
higher dimensions, but there have been several proposals for defining monotonic C func- 
tions by holography. In particular, in |2| positive-energy conditions in the bulk were 
used to produce a monotonic boundary function in terms of the superpotential of the flow. 
It is now interesting that in our two-dimensional example, we have both these objects at 
our disposal: Zamolodchikov's C function in terms of the holographic correlators, as well 
as the holographic proposal of 12] • 

The paper is organized as follows. In section 2 we review the exact domain wall 
solution of [5] that describes the flow. In section 3, we display the field equations and 
solve them perturbatively, i.e. close to the AdS boundary. In section 4, this solution is 
used to compute counterterms to form the renormalized action, which is then functionally 
differentiated to give the one-point functions. Solving the bulk fluctuation equations, this 
is already sufficient information to determine the two-point functions for inert scalars 
in section 5. In this section we also address the conceptual issues of prepotentials and 
distinction between A + and A_. For the sector of active scalar and stress-energy tensor 
fluctuations, one needs to do some more work, since their mutual coupling requires one 
to find gauge invariant quantities to work with. This is done in section 6 and we compute 
their linearized fluctuation equations around the domain wall. Then, in section 7, we use 
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the stress-energy correlators to study the Zamolodchikov C function and contrast it with 
other proposals for a C function. 

Since this paper is relatively long, let us give some guidelines on how the reader can 
get the most out of it in the shortest amount of time, depending on his or her preferences. 

• Learn the formalism of holographic renormalization: this paper provides an example 
that is technically less demanding than in the defining papers jTQl HH H2] , where 
emphasis was on the 5d/4d case, yet is still quite different from the basic example 
of a free massive scalar in the review ^3]. The reader with this interest would be 
well-advised to concentrate on sections El and EJ the philosophy is explained in the 
latter, with the core explanation in section l4~Tl Then it is straightforward to derive 
inert correlators as in 15 A\ the active / metric sector also requires decoupling their 
fluctuation equations, as in section El 

• See new conceptual issues that are also relevant in other dimensions, like the inclu- 
sion of inert scalars, the discussion of prepotentials, and the distinction between A + 
and A_ using correlators. These issues are dealt with in section El The notation is 
fairly standard in the literature so skipping previous details should not encumber 
the expert reader much. 

• Readers interested mainly in deformations of CFT, not supergravity details can 
enjoy our main results, the deformed 2-point functions, in section 15.41 and at the 
end of section 16.21 The notation is mostly self-explanatory, but it might be wise to 
simultaneously consult appendix B where ^ a (p) and the relevant special functions 
are explained. The C function is presented in section [71 

2 An exact holographic conformal-to-conformal flow 

In this section, we briefly review the analytic domain wall solution of [3] which interpolates 
between two AdS vacua. It was constructed as a solution in the three-dimensional N = 8 
gauged supergravity with local SO(4)xSO(4:) symmetry jB] , describing the AdS 3 xS 3 xS 3 x 
S l near-horizon geometry of the double D1-D5 system |2TH I7j with equal D5-brane charges. 
The matter sector of this theory consists of n multiplets each containing 8 scalars and 8 
fermions, whereas graviton, gravitini, and the 12 vector fields are non-propagating in three 
dimensions. The 8n scalars parametrize the coset manifold SO(8,n)/(SO(8)xSO(n)). 

The supergravity Lagrangian is given by 

C = ^VGR + Ccs + j^/GG^V^V^ + VGV + Cf, (2.1) 



5 



where Cqs is the Chern-Simons term for the vector fields, the third term is the kinetic 
term for the scalars, given explicitly in eq. ()3.3j) below, V denotes the scalar potential 
and finally Cp contains the fermionic terms, given in [5]. We use indices I, J, . . . and 
indices r,s, ... to label the vector representations of 50(8) and SO(n), respectively. The 
12 vector fields transform in the adjoint representation of the gauge group 

S0{4) + x S0{4)- c S0{8) c S0{8) x SO(n) C S0{8,n) , (2.2) 

where we use superscripts ± to distinguish the two three-spheres. In addition to the 
local gauge symmetry, the theory is invariant under the rigid action of SO(n). Assuming 
n > 4 matter multiplets, we break the latter down to SO (A) x SO(n— 4) and consider the 
following subgroup 

G inv = S0(4) inv x SO{n-4) 

C (S0(4) + x S0(4)~) x (S0(4) x SO{n-A)) , (2.3) 

of the global invariance group of the potential V. The S0(4) inv factor in G inv is embedded 
as the diagonal of the three SO (4) factors on the right hand side. Evaluation of the scalar 
potential V on the two-dimensional space of singlets under Gi nv leads to the potential jHj 

V = -g 2 (lQ + 24(Zl + Zl) + 8(Zl + Zl) 2 -8(Zl + Zl)-4(Z* + Z*) 2 ^ . (2.4) 

where g is the remaining gauge coupling constant. 




-1 -0.5 0.5 1 



Figure 1: Contour plot of the scalar potential V(Z\, Z%) ()2.4|) and the flow trajectory. 

The form of this potential is depicted in figure ^ It exhibits two inequivalent extremal 
points apart from the local maximum at the origin. The saddle point at (Z ly Z 2 ) = (1, 0) 
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corresponds to a nonsupersymmetric but stable vacuum as has been verified by explicit 
computation of the scalar fluctuations around this point [3]. We will concentrate on the 
extremum located at (Zi, Z 2 ) = (1, 1), which preserves N = (1, 1) supersymmetry. 

The ratio of the central charges of the dual conformal field theory at this extremum 
and that of the CFT at the origin is given by |2*T| 12*2*] 



cm 
cuv 



Vuv 



1 
2 



(2.5) 



supporting the conjecture that this point corresponds to a mass deformation of the UV 
conformal field theory; half the fields are integrated out to form the IR theory. The 
supergravity spectrum around this point is organized in N = (1, 1) supermultiplets as 
summarized in table HI where h and h denote the conformal dimensions associated with 
the supergravity masses, so that A = h + h. Note that the multiplet in the (1,1) for 
instance contains two scalars and spin-| fields, whereas the multiplet in the (1,3) combines 
a spin- 1 field with massive self dual vectors and a massive gravitino. 



SO(4) inv 


N — (1, 1) multiplets 

h x h 


field content 
(h,h) 


(1,1) 


(I 


i) x (III) 


( 5 5 ) ( 7 7 ) ( 5 7 ) ( 7 5 ) 

V4?4/? V4'4/' V 4 ' 4/? V4? 4/ 


(3,3) 


(i 


I) x (III) 


f 1 *) ( 3 3 ) f 1 3 ) ( 3 M 

V4?4/? V4'4/' V45 4/? V4? 4/ 


(1,3) 


(1 


i)x(fii) 


f 1 ? ) f 1 5 ) ( 3 7 ) ( 3 5 ) 
V454/5 V45 4/? V4' 4/ 


(3,1) 


(I 


i) x (ill) 


C -) ( 5 \) ( 7 3 ) ( 5 3 ) 

V4?4/? V 4 ' 4 > ' V45 4/? V45 4/ 


(2,2) 


(1 


1) x (||1) 


(1,1), (1,1), (|,1), (1,|) 



Table I: Supergravity spectrum around the supersymmetric vacuum (Z\, Z2) = (1, 1). 

An analytic domain wall solution interpolating between the origin and the supersym- 
metric extremum was constructed in [Sj using the ansatz 



Zi = Z 2 = -i=sinh(^g 



(2.6) 



where Q denotes the active scalar field parametrizing the diagonal in figure The 
remaining scalar fields are collectively referred to as inert scalars. Around the origin 
Z\ = Z 2 = 0, the mass of the active scalar field is m 2 L^ = — |, i.e. it is dual to a relevant 
operator of conformal dimension A = | which drives the flow away from the UV conformal 
field theory. In the truncation (|2.6)l . the bosonic part of the Lagrangian ()2.1j) reduces to 

£ = I VGR + I VGG^ d^Q d u Q + V Q , (2.7) 
with a scalar potential derived from a superpotential W as 

V Q = ^ {d Q W) 2 - 2W 2 , W = -| (l3 + 20 cosh(v / 2Q) - cosh(V8( 



• (2- 
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With the standard domain wall ansatz 



ds 2 = e 2A{r) rjij dx i dx j + dr 2 , 



(2.9) 



the field equations may be reduced to first-order form: 



d r Q 



d Q W, 



d r A 



- 2W . 



(2.10) 



They may analytically be solved by 



(5-y)(y+l) 2 
16(y-l) 3 



2igr 6A(r) 



(5 - y) A 



y = cosh(V2Q) .(2.11) 



128 (y+l)(y-l) 6 ' 



This solution interpolates between the origin and the supersymmetric extremum, preserv- 
ing N = (1,1) supersymmetry throughout the flow. From now on, we will fix the gauge 
coupling constant to the numerical value g = 1/8, thereby setting L , the AdS length at 
the origin of the scalar potential, to unity. 

Let us emphasize that for the computation of correlation functions, it is indispensable 
to have an exact, meaning analytic, domain wall solution (|2.11jl . rather than just finding 
the first-order equations and solving them numerically, which may be sufficient for other 
purposes. The point is that it is the fluctuations of supergravity fields around the domain 
wall solution that contain information on correlators in the boundary field theory. In 
particular, these second-order fluctuation equation are supplied with boundary conditions 
that specify the value of a supergravity field at the AdS boundary and demand a given 
(regular) behavior in the deep interior of the bulk space. The latter condition cannot 
be fixed perturbatively from the boundary, but is fixed only by solving the fluctuation 
equations. One might say this requires the "nonperturbative" bulk information that 
encodes the two-point function. 3 The fluctuation equations we obtain are highly singular 
at infinity, and it is not clear that one could have obtained our correlators even numerically 
if our domain wall solution would have amounted to just a numerical solution of the first- 
order equations. 

3 Field equations and near boundary analysis 

Thus, we embark on the way to computing correlators using the domain wall solution 
given above. As a first step we will derive the supergravity field equations and determine 
the coefficients in the near-boundary expansion of the supergravity fields. The truncated 

3 Of course, "perturbative" here refers to the radial expansion, which is a derivative (low-energy) 
expansion, rather than the coupling constant expansion in the field theory. 
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supergravity Lagrangian is given in (|2.7|) . In principle, it contains the complete informa- 
tion to compute correlation functions of the stress-energy tensor along the flow. Since 
we will be interested in computing correlation functions of the operators associated with 
the inert scalars as well, we extend (|2.7|) by expanding the original Lagrangian (J2.1|) to 
second order in the inert scalars. This is sufficient for all further computations in this 
paper, as we will only treat correlation functions with at most two insertions of inert 
scalar operators. We denote the inert scalars collectively by 

$ i = , (3.1) 

where the index i denotes the dimension of the representations under 50(4);^ while the 
superscripts ± label the two-fold degeneracies, cf. the spectrum in table |TJ Recall that 
the scalar fields in the (3, 1) + (1, 3) appear in a multiplet together with the vector fields 
to which they are related by gauge symmetry. Hence, these scalars require a separate 
analysis and will not be treated in this paper. 

The scalar fields parametrize an 50(8, n) matrix according to 4 

S = S Q St>, with 5$ = exp^$ i r i , (3.2) 

i 

where Sq carries the entire dependence on the active scalar Q, and Y 1 denote the non- 
compact generators of 5*0(8, n) associated with the representations of see (HI E] for 
details. With this ansatz, the current in the kinetic term of ()2.1|) becomes 

V lrylr = <j-l = ^^^1(^)5^ (3.3) 

where the /, r indices belong to the vector representations of 50(8) and 50 (n), respec- 
tively, and the relation to the i representation sector index is as described in section |5J 

The near-boundary analysis is most conveniently performed in Fefferman-Graham co- 
ordinates (x\p = e~ 2r ), taking the metric to be 

ds 2 = — - dp 2 H — g i jdx l dx :i i — 1, . . . , d , (3.4) 
Ap 1 p 

where for the moment we let d, the dimension of the boundary, be arbitrary for the sake 
of easy comparison with standard literature. The action on the asymptotically AdS space 
M is then of the general form 

S = [ d d xdpVG(±R+\G^K(<b)d^Qd v Q + \G^ y"d^%& + V(Q,$) 



4 To avoid confusion, we will always explicitly write out the sums over the different representations 
labelled by i. 
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d d x^lC 



(3.5) 



dM 



for general gravitational coupling constant k. We denote by K, the trace of the extrinsic 
curvature tensor tC^ on the hypersurface dM, and /C^- itself is given by 



Ki 



9ij dpQij 



(3.6) 



in the metric (|3.4|) . The induced metric on the hypersurface dM is denoted 7^. The 
functions V(Q, $) and i^($) in ()3.5|) are obtained from expanding (to second order in $) 
the scalar potential and the kinetic term, respectively, of ()2.1|) . 5 Somewhat miraculously, 
the quadratic parts of the two functions may simultaneously be diagonalized with Q- 
independent eigenvectors in each of the two-fold degenerate representation sector 4 and 9, 
respectively, and take the form 



V(Q, $) = V Q + J2 Vi(Q) , K($) + 



(3.7) 



with 



Ki = K 9 - = K^ = l 



Ka 



(3.1 



and the V\(Q) given in (|5.1jl . (|5.2|) below. The equations of motion from varying this 
action are, for a general metric G 



R/j,i/[G) + 



[IV 



— K 



2G„ 



Kd^QduQ + ^9^%^ + j^j (V - V (0) ) 



d 



<9$ 



x (V+\KG^d^Qd v Q) 



(3.9) 



where V( ) = V(Q, $)|q = o,$i=o, and A = «V(o) = — q!( - c '~ 1 - > , and @gQ is the expression 
obtained by varying \\/GKG^ u d^Qd u Q and discarding a boundary term. With the met- 
ric ([3 .4)) . the Einstein equations can be written in the following form, 



p[2^ - 2(g'g- 1 g% + Tr (fTY)^] + - (d-2) <4 - Tr (<rV)# 



5 It is due to our parametrization (|3.2|l . (|3.3() that the inert scalars automatically arise with a canonical 
kinetic term, whereas the kinetic term of the active scalar Q depends on the inert scalars. In contrast, 
the commonly used parametrization (which corresponds to choosing S = S$Sq in l|H.2|l) gives a canonical 
kinetic term for the active scalar and a Q-dependent metric for the inert scalars. This requires additional 
rescaling of the inert scalars in order to diagonalize their equations of motion, cf. |23| , which essentially 
amounts to reverting to the parameterization i|3.2|) . 



3 Our curvature conventions are R u 



8 V T 



1 1 ua 



(jj, «-» v), = Rfj 
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— K 



VtTr (g- l g')-V j g' 



Tr (<rV') - ±Tr OrY^ V) 



i 



-2k 
-2k 



(3.10) 



V-Vt 



(o) 



2(d-l)p 2 



and we note that here Rij [g] means the Ricci tensor for the <i-dimensional metric only, 
4pV" + 2p$ i/ [(2-d) + p (log (7)1 + pD^ 1 



i.e. only Rikj h and no -Rj P j P piece. The scalar equations take the form 



(2V, + | K^Q'Q' + pg^diQdjQ)) & , 

4p 2 (W + ^ , g / ) + 2pW[(2-rf)+p(log(?) , ] + p@ 9 Q-9QV = 0. (3.11) 

Primes here denote derivatives with respect to p, and g = detg^. From now on, we 
specialize to d = 2, but at times we will return to compare our results to those of higher 
dimensions, d — 4 in particular. For reference, it will be useful to have the first Einstein 
equation for = f(p) 77^: 



2p 2 f"-2pf + 2k(V- V ( o))/ = 



(3.12) 



and the scalar equation of motion for x-independent Q = Q(p) in this metric, with $' = 0: 



3.1 Perturbative solution 

We will now solve the above equations of motion perturbatively for the first few terms 
in the expansion in p. First, we note that the kink solution (j2.11j) admits an expansion 
according to 



Qb 



p l/4 + 1 /p_il p+ .. 

H \ 24 VF 64(T 



(3.14) 



where "B" is for background. Consistency of these expansions can be conveniently checked 
using (|3.12|) and (|3.13|) . The appearance of square roots of p in these expansions is generic 
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in three bulk dimensions. One could have chosen to expand in a different variable to 
avoid the square roots, but we find it more useful to stay notationally close to the higher- 
dimensional literature. In general, the metric and scalars may also depend on x, and we 
expand the fields in p as 

Gij(x,p) = p^g^x.p) (3.15) 
= P" 1 {g(o)ij(x) + y/P9(i)ij{x) + P (9(2)ij(x) + h {2 )ij{x) logp) + . . . 

Q(x,p) = p 1/4 g(x,p) = p 1/4 |g (0 )(x) + v /p(g {1) (a;)+g ( i ) (a;)logp) +...| , 

$(x,p) = p 1/4 0(x,p) = p 1/4 |0 (O )(x) + VP (^(l) (^) + 0(1) W logp) +...} • 

The subscripts in parenthesis denote the order in v /p, and also the highest number of 
derivatives with respect to x that will occur in these coefficients, hence this is a derivative 
(low-energy) expansion. The coefficients of the leading terms g(o)ij, Q(o), 0(o) are inter- 
preted as source terms for the dual operators. The field equations then determine the 
next few coefficients as algebraic functions of these boundary data. At a certain order 
- here, g^) for the metric and 0(i) for the scalars — the desired coefficient cancels out 
of the perturbative equations of motion, and remains undetermined; it is related to the 
one-point function of the dual operator. At this order, the ansatz is generalized to include 
a logarithmic term as we have done in ([3. 15)1 . and the coefficient of this logarithmic term 
is determined instead of the one that cancelled out. In general, for example in the d = 4 
Coulomb branch (CB) flow (see e.g. p]), one may also need higher powers of logarithms 
in the ansatz, but we have checked that we do not. 

The expansion of the scalar potential Vq gives 

13 1 3 

Vq = ~2 ~ g1{o)VP- g(^ )+ 6 ?(o)9(i))p-^g(o)g(i)plogp + --- , (3.16) 

while for the inert scalar potentials V\ we find 

Vi = -l-\(l + \K i )qf 0)v fp-+... . (3.17) 

To this order, the potentials V\ come in just two different forms, depending on whether 
the value of K\ (13.8)1 is or 1. It is only at higher order that the potentials V\ begin to 
differ between the various inert scalars. We can now expand the full potential as 

V(Q, $) = V {0) + 4p 7(1) (a:) + v/p logp V {1) (x) + p V {2) {x) + plogp V {2) (x) + ... 

and find for the coefficients 

V(o) = -\ , V m = - I (gf 0) + J^o)) 2 ) , V {1) = , 
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(3.18) 



i 



To lowest order, only the potential Vq for the active scalar contributes. The fact that Vn) 
vanishes is required for consistency; there is a 1/p in front of V in the equation of motion, 
giving a total of p _1,/2 logp for the Vm term, but there is no term of that order to match 
it on the left-hand side of Einstein's equation. 



3.2 Metric coefficients 

Solving the (ij) component of the Einstein field equations (|3.1Uj) for each coefficient in 
flSZEH) leads to 7 

9(i)ij = |«V(i) 9(o)ij , Tr# (2) = \R[g{Q)\ + 2kV {2) + |« 2 Vgj , 

h(2)ij = kV{2) g(o)ij , (3.19) 

at orders p -1 / 2 , p°, and logp, respectively. When one arrives at order p° of the Einstein 
equation, which would determine the coefficient g(2)ij, this coefficient only appears traced 
on the left-hand side. Thus, the non-trace part of g(2)ij remains undetermined in per- 
turbation theory, which is expected as remarked above. Note that the last equation in 
f!3.19|) is understood to hold only up to terms in order 0(L, since our starting action ()3.5|) 
was valid up to this order, cf. the discussion in the beginning of this section. It is further 
worth pointing out that, in contrast to the known d = 4 examples (e.g. the GPPZ flow), 
the h(2) coefficient in ()3.19j) is excited only by logarithmic terms in the scalars. 

From the (ip) component of Einstein's equations, one further derives 

3 1 
V J g m = Vi Tr g {2 ) - - V, Tr gf 1} - - g m V J Tr g {1) 

1 K 

+ 2 VJ + 2 feo) V* q m + 3 g (1) Vj q {0) ) 

+ \ E 0(0) V * ^(i) + 3 ^(i) V ^ ^(0) + | W( )) 2 <7(0) V, q {0) ) . (3.20) 

i 

This expression will be used in section (14 .5j) to verify one of the Ward identities. 

7 Here and in the following, we use the matrix notation tr A = tr (g^A) and AB = Ag^B, i.e. indices 
are raised and lowered with the boundary metric g(o)y- 
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3.3 Scalar coefficients 



The equations of motion for the inert scalar fields ()3.11|) in d = 2 read 

2p 2 &" + p 2 (\ogg)'&' + \ P U g & = (V i + K i (f?Q , Q , + \pg*a i Qd j Q))<!> 1 . (3.21) 
Expanding the right-hand side of (j3.21j) in p, we find 

- \p l " (#o) + + Mo)1(o))VP + 4>\i)VP^gp + • • •) • (3.22) 

Remarkably, K\ drops out of the expression, i.e. to this order in p, the effective potentials 
of all the inert scalar fields coincide. Now expanding also the left-hand side of (I3.21J) . we 
find that equating the most divergent terms (orders p -1 and p _1 / 2 logp) on the two sides 
just gives consistency conditions. The next-to- leading-order nonlogarithmic terms p -1 / 2 
have the same dependence on on both sides, so it cancels out. With the logarithmic 
term in the ansatz (j3.15|) . the coefficient 0m also appears at order p -1 ^ 2 , but it does not 
cancel and is determined as promised: 

<%) = i(Tr( 7(1) + 2g 2 o) )0 i (o) = 0{<f ) , (3.23) 

which upon using (J3.19j) is of only cubic order in the inert scalars, i.e. it vanishes to the 
order of validity of our computation for all the inert scalars. As anticipated above, the 
coefficient 4>\\) remains undetermined by the field equations. For the active scalar Q, one 
derives from ()3.11|) the analogous equation 

9(D = -i(Tr^ (1) + 2g ( 2 o) + 2^(0 i (o) ) 2 )g (O ) = 0(0 4 ) , (3.24) 

i 

while g(i) remains undetermined as expected. We are now in a position to use these results 
to compute divergences of the on-shell action and, from there, counterterms and one-point 
functions. From now on, we set n = 2 for the gravitational coupling constant. 

4 Counterterms and one-point functions 

4.1 Counterterms 

As is well known, the on-shell action of gravity and scalars is a priori divergent even on 
a background of pure AdS. It contains divergences coming from the p> — ^ (boundary) 
side, and the standard regularization is to consider the action on a surface away from 
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the boundary p = e for some small coordinate distance e. Then, divergences appear as 
poles and logarithms in e that can be subtracted, which corresponds to some subtraction 
scheme in the boundary field theory. Finally, one can let e — > to obtain the renormalized 
action. In addition to this standard procedure, we also want to ensure the on-shell action 
vanishes when evaluated on the background; this is necessary (but as we will see in section 
not sufficient in our case) to make sure the renormalization scheme is supersymmetric 
JBE|. Since this formalism is not yet widely familiar, we will make an effort to explain 
all steps carefully. (In this section, for transparency of the intermediate expressions, we 
suppress the indices i labeling the representation sectors of inert scalars, and restore them 
at the end.) To start from the beginning, the first step of holographic renormalization is 
one of convenience: to eliminate the Einstein-Hilbert term in the bulk action against all 
scalar kinetic terms on-shell. This is an amusing exercise that can be performed in any 
dimension, simply by tracing the Einstein equation in (J3.9)) and substituting the Ricci 
scalar so obtained into the action ()3.5j) . Only the potential V then remains in the bulk 
action, with the coefficient changed from +1 to —2, and there is a factor of 1/2 from 
y/G = v fg/2p 2 . Explicitly, 

S reg = / d d xd P VG (±R + \G>™K{$)d ll Qd u Q 
Jp>* v ^ 

+ L G ^J2d^%& + V(Q,$))-± / d d x^)C, (4.1) 

i ' J p= e 

= -j j P "d 2 xdp^r V -\j p d 2 x^-K, (4.2) 

where in the extrinsic curvature term we replaced yfy = y/g/e. Although expressions 
(|4.1j) and (|4.2|) are on-shell equivalent, the action ()4.2j) is only convenient for determining 
counterterms and later we return to using (|4.1j) . Here, the upper limit of integration p cr (x) 
can be defined in general as giving the surface p = p cr (x) of vanishing area, as in [24 , but 
for a stationary metric as we consider here, all one needs is that the bulk lapse and shift 
functions vanish at p — p CT . In the coordinates ([3.4)1 . this simply means p CT = oo. 

One then performs an expansion in e of ()4.2j) . Gravity without scalars in three bulk 
dimensions is expected to give a constant volume divergence plus a finite term involving 
Tr #(2). However, when we couple scalars we get all kinds of terms, including (logp) 2 due 
to $ 2 , as we shall see. We will use the expansion of y/det g to finite order: 

V9 = { 1 + v^(|Tr g (l) ) + e(|Tr g {2) ) + e log e (|Tr h {2) ) + O^' 2 ) } . (4.3) 

Notice that the combination Trg 2 ^ — ^(Tr c/(i)) 2 vanishes for our solution, since g(i)ij is 
simply proportional to g(o)ij- Substituting in all expansions, and staying with a general 
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potential, we arrive at 



S, 



reg 



j d 2 x ^ {e-\-V {0) - 1) + e-^i-Trg^V^ - 2V (1) - |Tr <? (1) ) 

+ loge(±Tr p (2) V ( o) + §Tr<7 (1) y (1) + % - ^Tr /i (2) ) 

+ (loge) 2 KV(2) + ^h m V [0) ) + e°(-^Tr^ (2) ) + ©(e 1 / 2 )} . 



(4.4) 



Of course we consider d = 2, but we found it useful to keep coefficients of the form (d — 2) 
to display some cancellations. In particular, we see that the finite Tt g^) contribution 
from the extrinsic curvature is cancelled. Substituting our perturbative solution for the 
metric (eq. (|3.19Jl ) into (|4.4j) . we stumble upon some pleasing simplifications: 



From a pure supergravity point of view, the cancellations in the loge term seem quite 
surprising: the back reaction of the metric to the scalars exactly cancel the explicit con- 
tributions of the scalars and leaves only the curvature. Also in the GPPZ flow [TT], there 
is no new "cross-term" anomaly, even though in that case as well as here, fluctuations 
of the active scalar and metric are coupled. In fact, holography provides a simple expla- 
nation for this. The boundary expectation value (T/), with scalar sources set to zero, is 
the trace of the variation of the generating functional in the boundary field theory with 
respect to to the boundary metric g(p)ij{ x ) — a l so with scalar sources set to zero. Hence, 
gravitational and scalar anomalies on the boundary can be separately computed, and 
must therefore simply add. On the other hand, one could have expected a separate scalar 
matter anomaly here like the $D 7 $ + |_R[7]$ 2 in four dimensions, but with R[y] = eR[g] 
and D 7 = eO g we see that these terms will vanish as e — > in our case. 

We now proceed to regularize the action with local covariant counterterms as outlined 
in [TTJ. The philosophy is the following: first, compute the divergent part of the regularized 
action S Teg for a given potential. Second, pull back all quantities to the regulating surface 
p = e, i.e. express the coefficients in terms of Q, and analogously for the inert scalars. 
The full fields Q(p,x) and x) are covariant, not the individual coefficients in the p 
expansion. Third, the covariant counterterm action S ct (Q,&) is defined by having the 
same divergent parts as S TCg (Q,<&) already obtained, and causing the final renormalized 
action 



S, 



reg 



/ 



d 2x V9(o){ e + e 



/ 2 (-f\/ (1) ) + log e (-| J R) + 0( e 1 / 2 )}. (4.5) 



S, 



ren 




(4.6) 



to vanish on the background, which implies (T^) = 0. 
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We now perform these steps in the present setting. For the inert scalars, substituting 
the V(i) given in ()3.18|) in the e _1//2 term, we find 

d 2 x^[ -- + + <j> 2 (0) ) - loge ±% (0) ] + 0{e 1 ' 2 )) . (4.7) 

To pull this back to p = e we need the inverse of (j4.3|) : 

^= v ^{l- v ^iTr^ 1) + e [i(Tr^ (1) ) 2 -|Tr^ (2) ]- e logeiTr/ i(2) + 0( e 3 / 2 )} .(4.8) 
Substituting flUJ) in flH7|) and again using (|3.19|) . we find 

Src g = J d 2 x4~g ( - 1 - i(g ( 2 ) + 0(0)) + (4-9) 

+e° (|(gfo) + 4)) - !^(o)0(D + f%(o)l) - loge |%(o)] + 0(e 1/2 )) - 

(Notice the sign switch of the term.) If we would add noncovariant counterterms 

to cancel only the divergences and then differentiate the remaining finite action, as one 
essentially did in the "old" prescription, we would find the wrong one-point function. 
What is wrong about it will be easier to see when we actually have the one-point function 
at hand, e.g. (|4.21|) for the active scalar. 

Instead, in holographic renormalization we concentrate on the divergences in the co- 
variant pulled-back action: 

S rcg = f d 2 x^ ( - 1 - i(I$ 2 + IQ 2 ) - loge \R\g] + O(e )) . (4.10) 

We remind the reader that it is 0(o) and g^o) which are kept fixed as e — > 0, hence e.g. 
Q 2 ~ e 1 / 2 as e — ► 0. Equation (|4.10|) determines the divergent parts of the covariant 
counterterm action S ct . In total, including an overall minus sign and using y/g~R[g] = 
the full 2d counterterm action is 

Set = /"rf 2 x v ^(| + iQ 2 + i^($ i ) 2 + aQ 4 + ^6 i g 2 ($ i ) 2 + loge| J R[7]) , 

"* i i 

(4.11) 

where the index i on $ has been restored. We did not only cancel the divergences in (|4.10j) . 
we also added the smallest set of finite counterterms needed for the renormalization scheme 
to preserve supersymmetry; the coefficients a and b\ are to be fixed. 

A few comments about the counterterms are in order. First, the most divergent term 
in the action is the volume divergence, which is cancelled by a counterterm of 1/2, or 
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(d — 1)/2L in general d and for arbitrary AdS scale L. This counterterm has been known 
since the early days of the AdS/CFT correspondence, for general d it was introduced in 
[23] . The quadratic counterterm |($ 1 ) 2 has also been familiar almost since the beginning. 
In |5] the coefficient of this term was computed in the fixed-background formalism for our 
case, and we indeed see it has the same coefficient (in general, (d — A)/2) as in the fixed- 
background case, justifying that treatment in retrospect. A finite Q A counterterm was 
first displayed in jTTj, and the Q 2 ($*) 2 terms can be seen as natural generalizations of the 
Q 4 term, with the additional complication that the coefficients b 1 cannot be determined 
by evaluation on the background, since $g = 0. This is not a problem; in section 15.31 
we will see how to compute b 1 , using a supersymmetry Ward identity for the two-point 
functions of superpartner inert scalars. The result is listed in ()5.27j) . 

Notice we could not add trilinear terms like Q($ 1 ) 2 since they would ruin the diver- 
gence structure; four powers of e 1//4 scalars are needed to reach finiteness. One could have 
added a finite counterterm proportional to the curvature scalar R, but as is well known, 
J d 2 x^fgR is a topological invariant in two dimensions (equal to Atc times the Euler num- 
ber) hence the variation of this term with respect to the metric #y is zero, so it would not 
contribute to correlators. 

Finally, we fix the coefficient a as in [TT] by evaluating S Ten = \im e ^ (S Teg + Set) on 
the background: 

S r0 n,B = J d 2 x (aqf Q)B + §tr (g(i) B )qfo)B + 5?(o)b?(i)b + Jtr (</( 2 jb)) 

= /A(a-±), (4.12) 

where we used the background expansions ()3.14j) . Hence a = 1/96, and we see that the 
naive subtraction of finite counterterms directly in the noncovariant action would have 
produced a different result. There is a useful check of the counterterm coefficients for the 
active scalar (see e.g. j2H]) : using the BPS equations, the expansion of the negative of 
W(Q) gives the coefficients directly, as -W(Q) = \ + \Q 2 + <^Q 4 + ... . 

4.2 One-point functions: generalities 

Now that we have the renormalized action, we can apply the formula at the heart of the 
AdS/CFT correspondence in the supergravity approximation: 




(4.13) 



field theory supergravity 



where the quantities on the left are the scaled Dirichlet data, e.g. q as opposed to Q. 
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The variation of the renormalized action for small variations in the sources g(o), 0(o) 
and c/(o) produces the boundary one-point functions: 

ss m = I d 2 x^ (U T <i) s 4) + W s m + EW^oi) • ( 4 - 14 ) 

Since bulk diffeomorphisms correspond to global symmetries of the boundary theory, one 
should require these boundary quantities to satisfy Ward identities. Following [TT] . one 
can derive Ward identities for the one-point functions of the stress-energy tensor and 
scalars in the presence of sources: 

v<r<,> = -{o q )v m - x;<e>;> v^oj , (4.15) 

i 

i 

= -| g(o) (O g }-|^0; o) (Oi) + ^, (4.16) 

i 

where A is the conformal anomaly we saw arise in the previous section, from a logarithmic 
counterterm breaking radial bulk diffeomorphisms. If the reader finds this form of Ward 
identities unfamiliar, it is probably because standard quantum field theory Ward identities 
are usually expressed with sources set to zero, so that e.g. V l (T^-) = 0. 

Explicitly, the anomaly A arises under a radial rescaling e — > /i 2 e; all terms in S ien 
above are manifestly invariant except the logarithmic term. It contributes an anomaly 8 
— ^\ogfi 2 A, hence we identify A = —\R. This is as was to be expected by holography; 
the conformal anomaly in 2d field theory on a space of scalar curvature R is simply 
proportional to R, 

L c 
-^ R= -2^ R = -2^ R={V) 

with c = 3L/2Gn the Brown- Henneaux central charge and the gravitational coupling 
temporarily restored to k — 8ttGn- The Ward identity (|4.16|) . including this anomaly, 
will provide a useful check for the one-point function of Tij computed below. 

4.3 Inert scalars 

The one-point function for the operators dual to inert scalars is now easy to compute, 
there is one contribution from the regularized action and some from the counterterms. It 

8 The factor -± is standard, cf. [TT] (5.43). 
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is convenient to also introduce the intermediate subtracted action S su b(e) = S vcg (e) + S ct (e), 
which is the quantity that becomes the renormalized action S iea when we take e to zero: 
S ren = lim e ^ 5 su b(e). The contribution from the regularized action (j4.1|) is 



5S reg = -2 I d 2 x^eJ2 S ® id e& 

Jp=e . 



where we used the $ bulk field equation ()3.9j) and y/g = Then, using the counterterm 
action S ct from eq. (j4.11|) we can write down the functional derivative 

' ' )KnU -- -2ed e & + \& + 2b i Q 2 & 



7 5& 1 

= eV4 [ - ¥U + ¥U + e 1/2 (-|0|i) + - 0V) + 2^(0)^0)) 
+ e 1/2 loge(-|0 i (1) + |0 i (1) ; 

This is divided by e A//2 = e 3 / 4 and the limit e — > is taken to yield the one-point function 
in the presence of sources 

(Op) = -(0 i (1) + 0| 1) )+26 i g 2 o) 0| o) = -0^ + 2^(0)0(0), (4.17) 

where we emphasize that this only holds to linear order in 1 , since we have only included 
$' up to quadratic order in the action. The forefactor of 0m is —1, or —(2 A — d) in 
general. (Notice that this is in ^T] conventions; there is an overall sign switch in the 
one-point function from JU]). This forefactor also reproduces the one obtained in j3] 
using the fixed-background formalism, where the finite counterterm that appears here 
was neglected. 

4.4 Active scalar 

For the active scalar, the contribution of the regularized action (j4.1j) is given by 

5S= I d 2 x^5Q(-2td e Q-e^K- 1 &&d e Q\, (4.18) 

Jp=e ^ ; 

where we used the bulk field equation for Q. Proceeding as above, we obtain for the 
one-point function 

(O q ) = -g (1) + 4agf 0) + g {0) ^(26 ; - \K,) . (4.19) 
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Expanding around the background (j3.14j) according to 



= «B(i) + ¥>(<) , (4-20) 

we see explicitly that vanishes on the background for a = 1/96 and there remain 
only the fluctuations: 

(O q ) = I^ (o) -^ (1) -I^ (o) ^(26 i -iirO0io)0?o) + ^ 2 )- (4-21) 

i 

This confirms that the flow is a true operator deformation; the boundary operator that 
drives the flow has vanishing vacuum expectation value. Notice that if we had not sub- 
tracted the finite counterterm aQ 4 in the action, or if we had tried to subtract a finite 
noncovariant counterterm aqf <. directly, we would have obtained a nonvanishing result 
for the one-point function (O q ); this would have been an apparent contradiction with the 
claim that the deformation has vanishing vev. 



4.5 Stress-energy tensor 

We proceed to compute the stress-energy tensor one-point function in the same way: 
functionally differentiate the renormalized action with respect to the induced metric 7^ 
on the surface p — e, but this computation is a little more complicated. It is convenient 
to split the computation into two parts, one corresponding to the extrinsic curvature, and 
one for the counterterms. The gravitational contribution due to the extrinsic curvature 
term is 9 



T 



rcg,grav, zj 



-\{-d e gij + gijg kl d e g M - 



p=€ 

d-1 



\ (JCij - Kj 



9ij) 



for general d, where by d e gij one intends d p gij(p,x)\ p=t . Here the factor of 2 in the first 
expression comes from ()4.14j) . but it immediately cancels with the 1/2 from the variation 
of ypy. For our metric ansatz we obtain 



K.'g.grav, 17 



9Wi) + e ~ 1/2 (~l9(i)ij + |#(o)yTrS(i)) 
+e (-2g {2 )ij - h {2)ij + ±g (1)ij Tr g {1) - ±g {0)lj Tr gf 1} + g {0)lj Ti g {2 )) 
+ loge{g m Trh {2) - 2h m ) + 0{e 1 / 2 ) 



'There is a sign mistake in (4-11). The correct expression is given in ^U] (3.6). 
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for the gravitational part. The counterterm action contributes 
2 5S ct 



X 



^7(57 



i.i 



~2 e 0(0)*/ 



I e 1/2 \_9m + 9(o)ij (9(0) + E^(o)' 



k(2)« + (?( 2 o) + E^(o)) 2 ) + l^(o)y (?(o)?(i) + E ^(0)^(1)) 



+ <7(o)ij- («9(o) + E 6i 9 (o) (4»> 

i 

/i(2)y + 0(D)ij (9(0)9(1) + E ^(0)^(1)) + ^(e 172 ) 



4 lo S e 



where the logarithmic term is scheme-dependent, being due to the variation of any matter 
conformal anomaly. Here, using the results of sect ions 13 . 21 and f3 . 31 all terms at order loge 
above actually vanish, but we will keep them since the expressions may be useful in 
situations where they do not vanish. Also, the Ricci scalar R does not contribute, as 
mentioned above. 

Putting everything together as T sub = T re g jgrav + T ct , the singular terms vanish upon 
using the earlier perturbative expressions. The finite part is the boundary one-point 
function, it is given by 



(Tij) = limT subi 

€— >o 



(4.22) 



§#(2)ij + \h(2)ij ~ j9(i)ij 1*0(1) + (9(b) + zJ^O 



o), 



+ 9(0) 



1.7 



i Tr #(i) " l Tr 9(2) - |9(i)9(o) - aqf 0) - ^ \ + 26; g 2 o) 0j o) ) <j>\ 



On the background this evaluates to 
1 







which vanishes by the value previously determined for a. This is, of course, a trivial 
consequence of the fact that a was defined to cause S Ten to vanish on the background. 

Using ()3.2()j) one may verify after some computation that ()4.22j) yields 



(4.23) 



VCfy) = (<Z ( i) - 4agf 0) - ^(26; - \K,) 0j o) 4)9(0)) V m 

i 

+E^(i)- 26i 9(0)^(0)) v^; 0) , 
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and thus satisfies the Ward identity ()4.15|) with the scalar one-point functions (|4.17j) . 
()4.19|) . Similarly, tracing (|4.22|) and using the perturbative expressions for Tr gm and 
Trp(2), we obtain 

(Ti) = -\R[g { o)\ + \ foi) - 4ag ( 3 0) ) q m + £ \ - (46, - 0(o) > 

i 

which together with f)4.17|) . (j4.19j) . guarantees the Ward identity ()4.16|) including the 
conformal anomaly A = —\R- Now we have our full collection of one-point functions and 
proceed to compute two-point functions. 



5 Two-point functions of inert scalars 



5.1 Fluctuation equations and the existence of prepotentials 

As has been emphasized above, near-boundary analysis is no longer sufficient when we 
move on to the computation of 2-point functions. It needs to be supplemented with a 
solution to the equations of motion linearized around the background (j2.11|) . For the 
inert scalars, it follows from (|3.11|) that these fluctuation equations turn into a three- 
dimensional Laplace equation in the domain wall metric, with total potential 

Vf* = V . + l Ki {d Q Wf . (5.1) 

Q=Qb 

These total potentials were computed in [5]: 



v^ ot = 


1 

1024 


(-45 - 160y+ 10y 2 + 3y 4 ) , 


T/tot 

^9+ — 


1 

16 


(17 + 30y + y 2 ) , 


T/tot 

V 9- — 


1 

1024 


(y + l)(-93 + 13y-19y 2 + 3y 3 


T/t0t 

V 4+ — 


1 

16 


(3 + y)(7 + 5y) , 


= 


1 

1024 


(y + l)(y-5)(l7 + 4y + 3y 2 ) , 



with y = cos1i(^/2Qb)- Further, it was shown in [Sj that the resulting equations of 
motion can be transformed into one-dimensional Laplace equations in flat space, with 
effective potentials Vi derived from prepotentials (in the sense of supersymmetric quantum 
mechanics) as Vi = U[ + U 2 . This underlying structure is crucial, as the absence of 
tachyonic fluctuations is then manifest. Although it is easy to see that V can always be 
rewritten in terms of a prepotential U for scalars with vanishing explicit potential V tot (i.e. 
when the effective potential V is only due to the e 2A of the curved background, cf. (J5.6j) 
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below), a priori it seems surprising that it would be possible for all our inert scalars with 
the various potentials above. A similar situation was observed in JS] for the active scalar 
fluctuations in the most prominent exact five- dimensional flows |HJ EJ 123, which seemed 
somewhat puzzling and in need of explanation. 

We now give a general argument for the existence of these prepotentials in the fluctua- 
tion equations of gauged supergravity and show how they may be directly extracted from 
the supergravity Lagrangian. Although we stay with our model as a concrete example, 
the argument straightforwardly translates to other supergravities and higher dimensions. 
The fluctuations of inert scalars around the background solution $' are described by the 
bosonic Lagrangian 

O = \^G%d^%& + y/G^V^(Q B ) , (5.3) 

obtained from ()3.5|) upon evaluation on the background. We change to horospheric coor- 
dinates 10 

ds 2 = e 2AB(c) fa,- dx i dx j + d( 2 ) , i.e. ^ = e~ A * . (5.4) 

dr 

Redefining $' = e~ A ^l 2 R} and dropping a total derivative, this Lagrangian takes the form 

C = \ d^R'd^R 1 + \ Vi R X K (5.5) 
in flat space, with a coordinate dependent potential 

Vi = 2e 2A »V?\Q B ) + K(C) + i(A' B (C)) 2 • (5.6) 
Restoring the fermionic part of ()2.1|) . cf. [H], we arrive at the Lagrangian 

C = I^W + iVi^ + lxy^ + l^ixV, (5.7) 
which is invariant under the global N = 1 supersymmetry transformations 

8R 1 = x l e , 6 X i = ^ed^ + UiR'e , 7 c e = e , (5.8) 

that explicitly descend from the Killing spinors of the domain wall solution (J2.ilj) . The 
fermionic mass term U\ now serves as a prepotential for the scalar potential 

Vi = U!+U? (5.9) 



10 Usually the horospheric coordinate is called z, we call it ( to distinguish it from the complex coordinate 
z we introduce later. 
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and may be extracted from the corresponding mass term in the gauged supergravity jH], 
more specifically from expanding the so-called A3 tensor around the background solution 
(j2.11j) to quadratic order of the inert scalars. The same prescription applies to all higher- 
dimensional supergravities. This shows that the existence of a prepotential is a direct 
consequence of the unbroken N = 1 supersymmetry of the background. 

Now, considering the full Lagrangian in the 9 sector, say, it is given by two copies of 
()5.7j) with prepotentials U\ related as IA 9 = —U 9+ . Closer inspection shows that on the 
total system we can realize another supersymmetry — i.e. in addition to ()5.8|) — with e of 
opposite chirality j^e = —e, and where the supermultiplets are (i? 9+ , X 9 ~)i (-R 9 ~> I n 
other words, the fact that not only N = 1 but N = (1, 1) supersymmetries are preserved 
implies that the potentials of the two scalars R 9+ , R 9 ~ are superpartners in the sense of 
super symmetric quantum mechanics, i.e. their prepotentials satisfy Wg_ = —Ug + . This in 
turn leads to the correspondence between solutions of the equations of motion 

R- = (d c -U + )R + , (5.10) 

which maps normalizable solutions into normalizable solutions. 

Having understood where the prepotential structure comes from, let us see how it can 
be exploited. To facilitate later contact to standard 2d CFT expressions, we switch to a 
complex coordinate z = (x 1 + ix 2 ) on the surfaces £ = constant. Then the plane wave 
ansatz is, using a complex variable also for the 2-momenta p = (p 1 + ip 2 ), 

R%z) =e i{p ~ z+pz) R\C) , (5.H) 
and the fluctuation equations take the form 

(-9 c 2 + V i ) J R i = -2|p| 2 i? i , (5.12) 
with the coordinate ( from (|5.4j) and effective potentials Vi derived from superpotentials [5j 
Wi = ±e A (y-l)(y+ll) , 

u Q ± = T^k eA (y -i)(y- 5 ) > 

U 4± = T ±e A (y-l)(y + 3) , (5.13) 

according to ()5.9|) . Curiously, the scalar <3> 9+ admits an alternative superpotential W g+ = 
—U-l, a circumstance that is not explained by the above argument; carrying different 
5 , 0(4) inv representations, $ 9+ and $ 1 can of course not be in the same (1, 1) supermulti- 
plet nor be related by another bulk symmetry. 
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Now, factoring out the zero mode (p = 0) solution, and switching to a new variable s 
related to y by y = 5 2 "Y~ 3 2 , the general solutions of the fluctuation equation ()5.12|) for R 9+ 
and R 4+ may be written as 

R 9+ = (i + y)- 1/4 x~2, 

with functions Xa(s) satisfying 

32 lol 2 

s X ^ + (l + a)x' a -^ L (2 + s 3 ) Xa = 0. (5.15) 

The remaining solutions are then automatically obtained by putting the supersymmetric 
quantum mechanics structure (j5.10j) to work: 

R 1 = (d C +Ui)R 9+ , 

R°- = (c\+W 9 _)i? 9+ , 

i? 4 ~ = (d ( +U 4 „)R 4+ . (5.16) 

The two ordinary differential equations (|5.15jl for a = 0, —2 thus comprise the entire 
dynamics of the inert scalar fluctuations. Remarkably, we will see below that the fluctua- 
tion equations for the active scalar and the metric reduce to the same universal equation 
f)5.15|l ! Moreover, even the fluctuations in the vector sector eventually lead to equations 
of type ()5.15j) [5 . This differential equation is a special case of the biconfluent Heun 
equation, and we analyze it and its solutions in appendix iBl 



5.2 Distinguishing A + and A 

The conformal dimension A of the operator dual to a scalar field of mass m in two 
dimensions is 

A = 1 ± Vl + m 2 . 

As first pointed out in P^j, there is a certain mass range for scalars (here —1 < m 2 < 0) 
where the negative root can make physical sense, hence one scalar can correspond to 
operators of two different physical conformal dimensions A±. Here we summarize some 
recent arguments why the distinction between the two roots is obvious for an active scalar, 
but is not for two inert scalars belonging to conjugate roots A± if they also have the same 
quantum numbers. 11 

11 Despite appearances, the "±" notation for e.g. 4+ and 4 ~ is not intended to imply that (/> 4+ has 
to be dual to an operator of dimension A+; in fact, we will see that </> 4+ is dual to an operator of 
dimension A_. 
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To begin, one can view the inert scalar expansions (j3.15|) as being composed of two 
independent interlocking Taylor series, beginning at orders p( d_A )/ 2 and p A//2 , respectively. 
The former can be thought of as the "source" (corresponding to standard AdS / CFT usage) 
series and the other the "response" series, in the language of [28J. In the series (13.15)) for 
A = A + , the exponent 1/4 is the (d — A + )/2, whereas the "response" Taylor series starts 
at order p A+ / 2 = p 3 / 4 and hence begins with the 0m and terms. As we have seen 
these terms are, in fact, independent of the sources 0( O ) and in near-boundary analysis 
(perturbation theory around p = 0), but acquire a "nonperturbative" interrelation upon 
demanding regularity of the solution in the bulk interior. 

For an operator with dimension given by the other root A_, in our case A_ = 1/2, 
the "response" series begins already at order A_/2 = 1/4, so essentially "source" and 
"response" are interchanged. Now, if one attempts to compute the 2-point function in 
the standard way, one is forced to add an additional quadratic counterterm, as explained 
in [2Ej- For a field with background (such as the active scalar q), this would generate terms 
linear in the source, which yields a nonvanishing 1-point function even when the source 
is set to zero. Hence, for the active scalar, selecting A_ is distinguishable from selecting 
A + in that the flows describe quite different physics, and the former choice presumably 
corresponds to turning on vevs in the boundary theory. 12 

For an inert scalar, however, these linear terms are not generated, and the only effect 
of the additional counterterm is that the coefficient of 0(i)0(o) in the action switches sign 
from — (A + — d/2) to(A + — d/2) = — (A_ — d/2). Thus, the flow and fluctuation equations 
are capable of describing the correlation functions of operators of dimension A_ equally 
well as A + . This means that it is not a priori obvious how to distinguish between inert 
scalars belonging to different roots A± if they also have the same quantum numbers, in 
our case 5*0(4) quantum numbers. In the next section, we will construct a supersymmetry 
Ward identity for two-point functions including finite counterterms, which will, among 
other things, allow us to make this distinction in section 15.41 

The 2-point functions for operators of dimensions A + and A_ are related by a "mas- 
sive" Legendre transformation, slightly generalized from ^7] to include a 0( O )^(o) term in 
the action with a coordinate- independent coefficient h^. As we have seen, such terms arise 
as finite counterterms, but unlike the coefficient a for the active scalar, h\ is not simply 
fixed by the domain wall solution (cf. the discussion after ()4.12|) ). In momentum space, 
the renormalized action for the inert scalar 1 can be written as 

SU4))) = § J 7050{o)(p)^o)(-P)(/+(bl) + 6l) , 

12 So far, however, the machinery of holographic renormalization has only been applied to flows whose 
active scalar is associated with an operator of dimension A + > d/2. 
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where /+ would have been the two-point function of operators of dimension A + , had we 
neglected the finite counterterm. The Legendre transform is effected by minimizing, just 
as in the massless case, the functional 



^(0(O)>0(1)) = £rcn(0( O) ) - / 7^2^(0) (p)^(I) (~P) 



with respect to 0'^, and in the linear approximation <f>\y = f + (\p\)cj)^ y One solves for 
the extremum in (f) 1 ,^ and substitutes in the original action to find 



In terms of the two-point function of the Legendre-transformed ("response") the- 

ory, this action should have the same sign as S(4>^), in analogy to the sign of the kinetic 
energy in transforming from Lagrangian to Hamiltonian for a free massive particle. We 
can thus identify the two-point function of the operator with dimension A_ as 

^ " " ' (5 ' 17) 

if is the would-be two-point function for the operator with dimension A + , when 

we neglect finite counterterms. This expression will be useful in section 15.41 



5.3 Counterterm fixed by Ward identity 

Now, what is this constant b\ and how can it be computed? To answer this, we will 
derive the supersymmetry Ward identity between correlation functions of inert scalars, 
embodying the supersymmetric quantum mechanics structure. In this subsection, let 
us again suppress the representation index i for clarity and only restore it at the end. 
Consider in general the correlation functions of operators dual to a superpartner pair of 
scalar fields ± . That is, the scalars defined by 

0± = C{p)B^ (5.18) 

for some function C(p) (in our case C(p) = e~ A ' 2 , but the explicit form is not important 
for the moment) satisfy the two fluctuation equations 

(-9 ( 2 + V ± ) J R ± = -2^1^, (5.19) 
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with V ± = ±(W ± ) / + (W ± ) 2 . The two (normalizable) solutions of (I5~TOJ) are then related 
by (j5.10j) . Let us first focus on R + , say. Expanding R + in a series in p, 



R+ = B+ +^B+ )+ pR+ +... , (5.20) 



1 f„l , ^ 



the correlation functions of the dual operators in the complex coordinate z = ^ (x + 
will be obtained in the next section as (see [T2*] ) 



(5.21) 



C(i) ^(i)^) 



C (0 ) R {0) (w) 

where C = C( ) + ^TpC^ + ... as usual. (Notice that in terms of the discussion of the 
previous section, this is for a "source" series; for a A_ operator, the correlator is actually 
the inverse of the right-hand side.) In the remaining part of this section we determine b. 

From d(/dp = (— 2pe A )~ 1 it follows that 

£ = -2^(1 + ^/2 + ...)* (5.22) 
<9C dp 

where the constant £ depends on the higher asymptotics in A(p), and is £ = — 1 in our 
case, but the precise value turns out not to be important. With the expansion ()5.20jh we 
find 

d ( R + = -^+-^2(2^++^+) + ... 

d 2 R + = 2R+ )+ £R+ ) + ... = R+ ) (V+ ) +2\p\ 2 ) + ... . (5.23) 

Here the last equality uses the fluctuation equation ()5.19|) and the expansion 

V = V ( o) + p 1/2 V ( i) + • . • =► V+ = - Ufa + (Ufa) 2 . (5.24) 

The latter equation follows from V = U + U 2 and the chain rule ()5.22|) . The expansion 
of the partner scalar R~ is then already determined from (|5.1()j) : 

R- = (d c -U + )R + 

= -J2+j - p 1 ' 2 R+ 0) (Vfa + 2\p\ 2 ) - (Ufa + p^Ufa^Rfa + P l7 Xi)) + • • • 

= + " P 1/2 {Ko) R t) + ( W (o)) 2 ^o) + 2 blXo)) + • • • ' 
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from which we obtain the ratio 



R 



(i) 



R 



(o) 



and using U 



(o) 



2 H 5 



R 

"<»> + fl 



(i) 



(5.25) 



(0). 



(0) 



U, Q y this gives rise to the supersymmetry Ward identity 




2 bp 



(5.26) 



for coefficients R^jRfQ) describing 13 two-point functions of two N = (1,1) superpartner 
scalars. Comparing to (|5.21|) . this determines the desired constant: b is simply given by 
^(o) — C(i)/C(o)5 i- e - b ± is given by the boundary value of the prepotential, up to a constant 
shift due to the curved background. Using the expansion ()3.14j) . one finds C(i)/C(o) = 1/4 
in this flow. We summarize the results as 



inert scalar 


U (0) 


h = Ul 0) - 1/4 


9+ 


1/8 


-1/8 


9- 


-1/8 


-3/8 


4+ 


-1/8 


-3/8 


4- 


1/8 


-1/8 



(5.27) 



These counterterm coefficients will be needed in the next section to show that the proper 
behavior of two-point functions emerges without dropping any terms by hand. 

As one would expect, the argument leading to the supersymmetry Ward identity (|5.26|) 
can be seen in a different way on a two-dimensional surface £ = constant, using known 
results in (1,1) superspace. Substituting d$R + = U + R + + R~ in (j5.7J) . we find that 
for one chirality choice of the parameter (rfe = — e), the supersymmetry transformation 
5x- has the superpartner R~ precisely in the place where the top component F of a 
(1,1) supermultiplet usually appears. In other words, the component transformations of a 
real scalar (1,1) superfield 5 with component fields (R + , x+, X-> R~) precisely reproduce 
(|5.8j) . Using standard results for relations between correlators of different components of 
the same superfield, one can finally recover a formula like (J5.26)) . 



5.4 Two-point correlation functions 

Since we study a flow to a fixed point, the boundary field theory is conformal in both 
IR and UV limits. This means we expect asymptotic power-law behavior of two-point 

13 Again, for the A_ operator the correlator itself is actually the inverse of the expression in parenthesis, 
as we will see explicitly in the next section. The Ward identity l|5.26|l then takes the familiar schematic 
form (OO) = 2\p\ 2 (00). 
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functions (the only exception being the 4 marginal scalars in the IR, cf. table [IJ) on 
both sides. In other words, (OaOa) — ► p 2A ~ 2 asymptotically, with two possibly different 
values Auv and Air. If we would apply the "leading nonanalytic term" (henceforth 
"old") prescription for AdS/CFT correlators, as successfully applied in e.g. [TH], we would 
quickly be disappointed in this case. This is because the correlator we are interested in is 
asymptotically just a power of p, and in the "old" prescription it would have seemed that 
the correct power-law behavior would be indistinguishable among other monomial terms, 
coming from unphysical contact terms, that were summarily dropped. 

In fact, in holographic renormalization, the requirement that only local counterterms 
may be added is sufficiently restrictive to single out exactly the right behavior. Indeed, 
yfyQy/U $ would yield a counterterm ~ p but is nonlocal, whereas V /7$D"$ would yield 
a p n counterterm but vanishes in the limit e — > for all our $ ~ e 1 / 4 scalars. In other 
words, holographic renormalization does not allow us to "drop contact terms" as in the 
"old" prescription, which is fortunate since precisely some of those monomial terms are 
physical ones in our context. (Of course, there are many examples where that prescription 
is still applicable and useful.) 

Now we compute the two-point functions of inert scalars. These two-point functions 
descend from uncoupled fluctuations that do not require an analysis of the kind we will 
perform for the active scalar and metric in the next section. They are obtained directly by 
taking the functional derivative of the one-point function (j4.17j) . and are thus essentially 
(i.e. up to finite counterterms) encoded in the ratio <fi}i\/<f>) \ in the expansion 

& = P 1/A (4>) + vWii) + • • •) , (5-28) 

of the solution to the fluctuation equations, after we impose regularity in the interior of 
the bulk. To compute this ratio, we recall that ^(Cz) = e i(p * +pz "> e~ A ^ /2 R l {C) , and 
-R^C) is given explicitly in terms of biconfluent Heun functions in (j5.14j) . (j5.16|) . Using 
the following asymptotic expansions of the domain wall solution (|2.9j) . (j2.11j) 

e -A/2 = p l/4 + I p 3/4 + A p 5/4 + _ ? 



y 



4' 64' 

l + p i/2 + I p _ JL p 2 + 
H V 192^ 



*-*o = \p 1/2 + r Q p+ v Q p3/2 + --- ■ (5 - 29) 

(where we recall that s is the variable used in the fluctuation equation (|5.15jl ). we even- 
tually find 

1 , 8\p\ 2 



8 + V- 2 (p) ~ 2 
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>0) 1 *- 2 (p) <i) 3 8|p| 2 



8 *_ 2 (p) 



3 tt (p) 00). 1 8b| 2 

<o) " 8 4 4fc " 8 *o(p) 

with the ratios ^> a from the expansion of the regular Heun function in (jB.2|) . Consulting 
table (|5.27jl . we see that the effect of the Q 2 $ 2 counterterms with the coefficients b\ 
determined above is precisely to cancel the additive constants in these expressions. The 
resulting correlators may be expressed in terms of the Heun function coefficients ^ a {p) 
(discussed in appendix IB]) : 

<0J(-p)0J(p)> = k *-'<">- ^ 



8(2 + 8|p| 2 - *_ 2 (p)) 



wrnow = (ot-(-P)ot-(p)) = 

and all mixed two-point functions vanish. Here we have restored the factor k/8 which 
comes from properly normalizing the supergravity action in (J4.13)) . This k is the level of 
the current algebra in the ultraviolet CFT, and it depends on the number iV of D5-branes 
as k ~ N 2 jjj. We only restore this factor in final results for correlators. 

The Heun function coefficients ^ a (p) are plotted in figures 3 and 4 in appendix iBl 
analytic expressions for the large and small p asymptotics are obtained in ()B.10|) - (jB.12|) . 
from which one may immediately derive the UV and IR asymptotics of these correlators. 
Appealing to (|5.17|) . we had to invert some of the ratios (|5.30p . namely those corresponding 
to the operators of UV conformal dimension A_ = 1/2. In the 4 and the 9 sector, we 
want to make the distinction between A + and A_. As it turns out, this distinction is 
unique for the scalars in the 4: it is only upon inverting the shifted ratio (<fify/<f>M\ — 3/8) 
that the resulting correlator has the correct logarithmic behavior from ([B.lljh as expected 
for a Air = 1 operator. Tracing this identification back through the flow, we conclude 
that in the UV it is <3> 4+ which is associated with A_ = 1/2, whereas $ 4+ corresponds 
to A + = 3/2. Thus, the computation of correlators throughout the flow allows us to 
make this distinction which would have been impossible to derive from a near boundary 
analysis around the UV boundary. 

In the 9 sector, on the other hand, the ambiguity is left unresolved. We have chosen 
to associate $ 9+ with A_ = 1/2 and $ 9 ~ with A + = 3/2 but could also have done it the 
other way round, both choices are compatible with the correct asymptotics of A_ = 1/2, 
A + = 3/2 operators in the UV and IR. 
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For the singlet $ 1 , we know that it corresponds to a A_ = 1/2 operator in the UV, 
since it is the superpartner of the active scalar which has dimension A + = 3/2 (see 
discussion in section I5~2")) . For this inert singlet scalar, we have fixed the coefficient of the 
Q 2 Q 2 counterterm such that the resulting correlator in (J5.31j) satisfies the supersymmetry 
Ward identity with the correlation function derived for the active scalar in (|fi.2()|) below. 
(It is quite a nontrivial consistency check that this Ward identity may indeed be satisfied 
by just adding the proper constant to the ratio in ()5.30|) .) This correlator leaves a minor 
puzzle in the IR: \I/_ 2 (p) goes linear in p for small p, so the correlation function goes to a 
constant rather than like the |p| 3 which one would have expected for A IR = 5/2. However, 
irrelevant operators decouple in the infrared, so it is difficult to know whether the method 
works straightforwardly for operators crossing over from relevant to irrelevant along the 
flow; certainly UV-irrelevant operators cannot be treated the same way as UV-relevant 
operators [TO] . 



6 Two-point functions of active scalar and stress- 
energy tensor 

6.1 Fluctuation equations 

Unlike the inert scalars, the active scalar fluctuation couples to the metric fluctuation at 
linear order in the equation of motion. This was a tough obstacle (see e.g. ^E]) until 
resolved in ^1] (with earlier progress in [29 .). The resolution involved working with 
"gauge invariant quantities", therefore we vow to only work with such quantities, in a 
sense to be made precise below. 

To compute the quadratic fluctuations of the active scalar and the metric, the inert 
scalars may be switched off as they do not contribute to these couplings. We use com- 
plex coordinates z = ^= (x 1 + ix 2 ) on the boundary, and parametrize the domain wall 
background and the fluctuations as 

ds 2 = e 2A{r) (2(l + h+\p\ 2 H)dzdz + p 2 (H+H ± )dzdz + p 2 {H-H ± )dzdz^ 
+ (1 + h rr ) dr 2 

q = q B + P (6.1) 
with the same plane wave ansatz in the complex z coordinate as was used earlier: 

(p(z, z, r) = e i{p ~ z+pz) <p(r) , etc. (6.2) 
for all the fluctuations. Linearizing the equations of motion around an arbitrary domain 
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wall (|2.9jl . we obtain 



^1 = 

i?"-4tra' = -e~ 2A h r 

2h rr W + ti = -AyW' , 

\p\ 2 e- 2A h = 2\p\ 2 W H' + (2<fJ -h rr W -2<pW")W , (6.3) 

in terms of the superpotential W(Q). In the above equations, primes on the superpotential 
W denote derivatives with respect to Q, whereas all other primes refer to derivatives with 
respect to the radial variable r. The first equation in (J6.3|) is a manifestation of the well- 
known fact that in three spacetime dimensions there are no transverse-traceless degrees 
of freedom in the metric. As explained in [IT], the ansatz (j6.1|) does not completely fix 
the bulk diffeomorphisms, but leaves the freedom of "gauge transformations" generated 
by vector fields (£ z ,£ 2 ,£ r ) satisfying 

d r C = -pe~ 2A C , d r C = ~pe- 2A C , (6.4) 

under which the fluctuations transform as 

8 h = 2f A' , 5 h rr = 2d r C , 8 <p = f V B , 
5H = ^ (pf- +P e*) , SH' = -2e~ 2A C, 

8H ± = -^(pC-pC), 5H' ± = 0. (6.5) 

It therefore seems appropriate to cast the above fluctuation equations in equations for the 
"gauge invariant" objects 

W Ae~ 2A 
Ji = h + A—<f } J 2 = -2H'-—<p, 

2(<f/ - cpW") 

K = h„ — . (6.6) 

This leads to 

M 2 e~ 2A Ji = WW - (AW 2 + {W'f - 2WW") K , 
\p\ 2 j 2 = (-2W" + AW) 1Z — 1Z', 

J[ = -2WK, (6.7) 

which we may combine into the second order equation for 1Z 

= K" + (2W"-8W)TZ' 

+ (2\p\ 2 e- 2A + 16W 2 -A(W') 2 -8WW" + 2W'W / ") TZ . (6.8) 
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The analogous equation in five dimensions was obtained in jllj . Specializing this equation 
to the superpotential W from ()2.8|) and transforming to y — cosh(y / 2<5), we find 



2\p\ 2 e- 2A TZ(y) = £(y - 5) 2 (y - lf{y + lfTZ"{y) 

+ Uv - 5 )(y- Viy + + % 2 )^'(?/) 

+ i (35 + y(80 + lOy + 3y 3 )) ftfo) . (6.9) 
Dividing out its zero mode as in (j5.14j) . 

K = (y-l)(y-5)-Hy + l)-$x, (6-10) 

and performing the same change of variables y = as that leading to f)5. 15j) . eq. ()6.9|) 
reduces to 



„ 4 + 5s 6 , 32 pi 

s * ~ YT^ * * (2 + ^ * = °' (6-11) 



We now proceed to show that also this equation may be reduced to one of the two Heun 
equations (|5.15jl . Indeed, let x be the regular solution of (j5.15|) with a = —2, then 

X = X-\s X ' (6.12) 

is a solution of (|6.11j) . i.e. we can extract all the asymptotics data from our previously 
obtained solution! In particular, after some computation it is seen that the solution of 
()6.8|) regular in the bulk interior has the expansion 

with \l/_ 2 (p) from f)B.2j) . Using this expansion we can now proceed to compute two-point 
functions. 



6.2 Two-point correlation functions 

The rest of the computation is straightforward. Linearizing the one-point function of the 
stress-energy tensor ()4.22j) around the background using 

9ij = fx? AB (rii3+hij) , ( 6 - 14 ) 
leads to 

( T ij) = Vij QTr h {1) - ±Tr h {2) + ^(o) - §¥>(i)) - + \h m , (6.15) 
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or, in the complex notation from (j6.1j) . 



(T,z) = \{h { x) + \p\ 2 H {1) ) - \{h {2) + \ P \ 2 H {2) ) + ^ (0) - , 
(T zz ) = f {-\H m + \H {2) ) , 

(T„) = p 2 {-\H {1) + \H {2) ) . (6.16) 

where the subscripts in parentheses denote coefficients in the p expansion, H = + 
y/pHm + pH(2) + . . . as usual, and we have used that H± = constant. Expanding the 
gauge invariant quantities ()6.6|) yields 



= (|^(o) - 4^(i)) 4f> + (-§ ^(o) + - 8^(2)) P + ■ ■ ■ , 

J{i) = (h( ) + 4^(o)) + + |^(o) + 4^(i)) y/p 

+ (h {2) + § y? ( o) + | ^(i) + 4^ (2) ) p + . . . , 

J m = 2H {1)y /p + (AH m + 8^ (0) ) p + • • • • (6.17) 

One sees that the coefficient lZ(o) = 1^(0) — ^f(i) is proportional to the one-point function 
(O q ) of the active scalar (just set (ft 1 = in (|4.21|) ). hence that the one-point function is 
gauge invariant as one would have hoped. From the equations of motion f)6 . 7j) . we then 
find that we can express all the perturbative coefficients in terms of v?(o)> ^(o)> an d the 
ratio 1Z(i)/1Z(q) which has been determined in (J6.13J) above: 

H n n i h , h {0) +4<p {0) 

— U , H(2) — g ft (0) + 



/i(l) = -2 V9( ) , h(2) = - \ </?(0) + 



4(ft (1) /ft (0) )-5 ' 

2\p\ 2 (h {0) +Aif (0) ) 



4(ft (1) /ft (0) )-5 

2b| 2 (fr(o) + 4y ( o)) 
A{K {l) /K m ) - 5 



Substituting this into (jOU), (|6~T5T) . we find 

2H 2 (/ i(0 )+4 v ? ( o)) 



4(K {1) /K {0) ) - 5 



1 <^-> = i<T„> = 1<IW = ^ (0)+ o/ >: 4 ?' 1n , (6.19) 



|p| 2 x 7 p 2 ' p 2 1 4 w 8(^(1) /K (0) ) - 10 
which gives the two-point correlation functions 

k\p\ 2 k\p\ 2 (-2 + ^_ 2 (p)) 



(O g O q ) 



A(n (1) /n m )-5 4(2 + 8|p| 2 -v]/_ 2 (p)) 
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(O q T zg ) 



-k\p\ 2 



k\ P \ 2 



-1 + \^- 2 {v) 



(6.20) 



WR {1) /K m - 20 



8(2 + 8|p| 2 - ^_ 2 (p)) ' 



k 4K {1) /K {0) -3 
32 lPl 4K {1) /K (0) -5 



k\ P \ 4 f 1 , 1 



8 V 8 N 2 2 + 8|p| 2 -*_ 2 (p) 



where we have substituted the expansion (j6.13j) for 1Z(i)/TZ(o). The first two correlators 
are, now that the smoke has cleared, trivially related by 14 



where (3 = —1/2 is the classical (3 function due to the classical scaling of the coefficient q^, 
when this coefficient is viewed as the holographic coupling constant in the deformation 
£cft + 9(o) O q . As in [TT], the fact that the (3 function is classical could be ascribed to 
a nonrenormalization theorem to the effect that the only contribution to the (3 function 
could come from an anomalous dimension of the operator it multiplies, and if this operator 
is protected, there are no quantum corrections to scaling. After all, the deformation 
preserves some supersymmetry, so this is perhaps not so surprising. 

It is straightforward to see that the two-point correlator of the active scalar has the 
correct (linear \p\) UV behavior, corresponding to an operator of conformal dimension 
3/2. Moreover, comparing to (|5.31)1 . we see that it indeed satisfies the supersymmetry 
Ward identity (O q O q ) = 2\p\ 2 (O^O^) . The third correlator, (T ZZ T ZZ ), and its asymptotics 
will be studied in more detail in the next section, since it is related to the C function 
along the renormalization group flow. 



Given the stress-energy 2-point functions, we may follow Zamolodchikov's original con- 
struction jTHj to compute the C function, a function on the space of couplings that is 
monotonic along the flow and interpolates between the central charges of the conformal 
fixed points. The variety of possible tensor structures in d > 2 makes this construction 
difficult to generalize to higher dimensions [3UJ l3"T| : in particular, the straightforward 
proof of monotonicity as a consequence of unitarity is tailored to fit the two-dimensional 
case. Nevertheless, there have been several proposals for defining monotonic C func- 
tions by holography [321 EH 12]- in particular the "holographic C function" of fHl |2] 
is a very simple proposal in terms of the supergravity superpotential W, which reads 
Choi = — ^/W(Qb) in two dimensions. It is monotonic as a function of the bulk radial 
coordinate, assuming a fairly weak positive-energy condition in the bulk supergravity. 

14 see e.g. [T3] eq. (6.13), and use (AO q ) B = 0. 



T zz = (30 q + A 



(6.21) 



7 The C function 
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Applied to our superpotential ()2.8|) . we find 
3k 



hoi 



W(Q B ) = S*(l-^ + ip + ...). 

with normalization adapted to our conventions. Monotonicity may be verified directly. 

We will now compute Zamolodchikov's C function from the holographic (TT) corre- 
lators obtained above. As a first check of the (TT) correlators ([6.20)1 . we expand their 
asymptotic behavior using (jB.10|) . (jB.12|) . and find 

<T„(p)r„(-p)> i -| ( i |^ + 2|p| J -i^|p| 5 + ...). (7.2) 

In particular, this shows that in fact cir/c^v = 1/2, which is precisely what is expected 
(see eq. (12.5)1 ). To construct Zamolodchikov's C function, we first Fourier transform the 
(T ZZ T ZZ ) correlator back to the complex z plane as 



(T z - z {z)T z -M) = i- J dpdpe^~ z+ ^ (T zz (p)T z - z (-p)) 



= : □ □ n(t) = (n" - n"' + in"") , (7.3) 

defining the function Q(t) with t = | log(/x 2 |z| 2 ), that is, <9 4 = |^|5| z |, and primes denote 
derivatives with respect to t. This function encodes Zamolodchikov's C function jTH| as 

^zam = -96 (n' - n" + Jn'") , (7.4) 

such that C = —24 \z\ 4 (T zz (x)T zz (0)) . Specifically we find the integral representations 

m = Jd\p\\p\- 3 J {2\pz\) (T z - z (p)T z - z (-p)) , 

n\t) = - J d\p\2\p\~ 2 J x {2\pz\) (T z - z (p)T z - z (-p)) , 

Czam = 192 J d\p\{\p\- 2 \z\ (HJ (2H) + (|pz| 2 - l)Ji(2|pz|)) 

x (T zz {p)T zz (-p)} } . (7.5) 

With proper regularization, the last integral gives 

db|{b|- 2 |^|(|^|Jo(2|^|)+ (7.6) 
(\ P z\ 2 - l)Jx(2\pz\)) | P r" 



irn\z 



2-n 



4sin(^f)r(2-f)r(-f) 
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for the polynomial terms in (T Z zT zz -). From the exact asymptotics ([7.2)1 . we may then 
derive the small distance behavior of the C function 

C Zam = 3k(l--± ? =\z\+0(\z\ 3 )) . (7.7) 



For the full C function we have to insert in ()7.5|) the complete expression from (|6.20j) . 
with \l/_2 from appendix B. The numerical result is plotted in figure El and it is another 
consistency check of the correlators (|(12())1 that the function comes out strictly monotonic, 
i.e. the holographic correlators indeed reproduce the properties of a unitary field theory. 

To compare the result to the proposal (|7.1|) above, we recall that in holographic flows 
the RG scale p is introduced by the AdS isometry 

p — > p 2 p , x % — > px % , (7.8) 



connecting the energy scale to the bulk radial coordinate. At fixed p this converts the 
superpotential W(Q&(p = p 2 po)) into a function of p. Likewise, Cz am turns into a function 
of p upon relating the boundary radius \z\ to a fixed value \z\ = p\z \. Normalizing Choi 
and Czam to unity in the UV and fixing the values of \zo\ and po such that the first 
derivatives C' hol and C Zam coincide in the UV, we plot the two functions in figure El 

i 

. 9 
. 8 
0.7 
. 6 



2 4 6 8 10 12 




Figure 2: Holographic C function ()7.1|) (solid) and Zamolodchikov's C function ()7.4|) (dashed) 
from holographic correlators as functions of p after normalizing and rescaling such as to match 
the first derivative in the UV. 



Both functions are strictly monotonic with surprisingly similar shapes. We should 
stress that the discrepancy in figure El does not mean that the different proposals are 
incompatible — the computation of the C function is in any case scheme dependent - 
but may rather indicate that the identification of the boundary energy scale with the 
radial AdS variable (J7.8)) requires corrections away from criticality. 
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A comparison of the different C function proposals in the five-dimensional confining 
flow of [1] was done in [33], see also [3^. Some general considerations in higher-dimensional 
conformal-to-conformal flows are presented in [35] . Note that in higher dimensions the rel- 
evant central function descends from the two-point correlators of the transverse-traceless 
part of the stress-energy tensor, which drastically simplifies the computation. In two 
dimensions one instead has to go through the full procedure of decoupling the active 
scalar and metric fluctuations, leading to the fluctuation equation ()6.8|) and its solution 
presented in the previous section. 

A challenging goal would of course be the computation of the (TT) correlators directly 
in the CFT as input to Cz am - Comparison to the correlators obtained by the holographic 
methods presented here could provide a demanding test of the correspondence. 

8 Outlook 

In this paper, we have computed correlators along a renormalization group flow interpo- 
lating between conformal field theories. As a main result, we have obtained the two-point 
correlation functions (TT), (TO), and (OO) of the stress-energy tensor, and operators 
O dual to supergravity scalars, both active and inert. We used the (TT) correlators to 
compute the Zamolodchikov C function and compare it to the holographic proposal for 
the C function in terms of the supergravity superpotential. Since we summarized the 
results in the introduction, here we only give some brief remarks on future directions. 

First, one could study correlators of currents J 1 in the boundary theory. As in |12j . 
these currents correspond to symmetries broken by the deformation q^o)O q (here SO (A)), 
although their analysis was facilitated by the remaining U(1)r symmetry of the IR theory, 
whereas our IR theory has no i?-symmetry at all. Presumably this analysis is nevertheless 
straightforward on the supergravity side, given the fact that also the vector fluctuation 
equations reduce to the same biconfluent Heun equation [S]. 

More interestingly, it would be worth to study the CFT side in more detail than has 
been done here. In fact, it should be possible to compute correlators such as (TT) given 
the deformation £cft + Q{o)O q , where g( ) is the holographic coupling. One could hope to 
show that this deformation is integrable, and to compute correlators exactly. 

One can also proceed to apply these methods in other settings. Perhaps the example of 
greatest direct physical interest would be domain wall solutions in five- dimensional gauged 
supergravity. Even though the 5d flow equations in j2] could not be solved exactly, one 
could see how far one can take a numerical analysis of correlators there. Unfortunately, as 
emphasized in section [21 it will probably be well-nigh impossible to achieve any confidence 
in the numerics if one cannot find any analytic help. On the other hand, it might be 
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possible to construct exactly soluble setups in other five- dimensional cases than that 
of [2 . Encouraged by the solution in section |21 and recalling that it was crucial for the 
explicit construction of [5 a to set the two D5-brane charges equal (other values of this 
ratio would have resulted in a vertical stretching of the potential in figure Q such that 
the flow trajectory would no longer be a straight line), one might look for simplifying 
special values of coupling constants in various five- dimensional gauged supergravities. 
In an analogous half-maximal theory in 5d, in which first-order equations for flows were 
studied in jSH], there is an SU (2) L x SU (2) R non-Abelian gauge symmetry but apparently 
no free parameter of this kind. 

As a matter of principle, it would be interesting to compute n-point functions for 
n > 2. Although the full Dirichlet problem may be practically unsoluble, one can imagine 
solving nonlinear fluctuation equations in perturbation theory. Recently, an outline of an 
example for n = 4 was given in [T3] . section 5.9. 

Finally, one could compute correlators in other RG flows in two-dimensional conformal 
field theories, further pursuing the idea of three-dimensional gauged supergravity as a tool 
in this field. Of particular interest and in principle accessible with our tools are for example 
flows to nonsupersymmetric but stable fixpoints, the analysis of marginal deformations 
of the CFT describing the D1-D5 system (cf. |3*T1 EH EH]), and flows in the maximally 
supersymmetric theory of [Hj and their role in a supergravity description of matrix string 
theory (cf. gDj). 

Acknowledgements 

We gratefully acknowledge helpful discussions with B. Farid, I. Runkel and K. Skenderis, 
and we especially thank M. Bianchi and M. Haack for useful discussion and comments 
on early drafts of this paper. We thank IHES, Bures-sur-Yvette, for hospitality when 
this work was being finalized. M.B. is supported by a Marie Curie Fellowship, contract 
number HPMF-CT-2001-01311. This work was supported in part by INFN, by the EC 
contract HPRN-CT-2000-00122, by the EC contract HPRN-CT-2000-00148, by the IN- 
TAS contract 99-0-590 and by the MURST-COFIN contract 2001-025492. 

A Notation 

Since we made many small notational changes with respect to [5], where the domain wall 
solution was constructed and discussed, it might be useful to collect the differences. The 
reason for these differences is the discrepancy between the 3d supergravity literature and 
holography literature (such as [H1H2]) an< ^ m this paper we stick close to the latter. The 
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differences between old (reference jS]) and new (this paper) notation are summarized in 
table El 



old 


Q 


q 


; 




Vi) 


9 2 V 




K 2 


new 


V2Q 




&/y/2 




(Zi,Z 2 ) 


V 


TAtot 

* i 


K 



Table II: Notation: reference [H] vs. this paper. 



In addition, the two integration constants in the domain wall solution (j2.11|) are 
rescaled relative to [S], so that also the expansions ()3.14|) appear different; in both cases 
the integration constants can be fixed by demanding g( ) B = 1 and g( )B = 1 in the domain 

wall solution. Finally, the signature in |H] is Lorentzian (H ) except for in section 4 

where it is Wick rotated to Riemannian signature, and here it is Riemannian throughout. 



B Analytics of the fluctuation equations 

We have seen in the main text that the entire set of fluctuation equations around the 
background (j2.11|) may be reduced to equations of the type 

s X l + (l + a)x' a -P 2 (2 + s 3 ) Xa = 0, (B.l) 

for a = 0, —2, and with P = a/32/3 \p\. This appendix is devoted to a closer study of 
this differential equation and the properties of its solutions. Let us recall that along the 
flow the variable s runs from s = 1 at the AdS boundary (the UV) to s = oo in the AdS 
interior. This implies that the two-point correlation functions are encoded in the first 
coefficient of the expansion at s = 1 

Xa(s) = + *«+•••) , (B.2) 

of the solution Xa regular in the interior s — ► oo. The coefficient is uniquely defined 
as a function of the parameter P by the requirement that Xa is regular as s — ► oo. This 
procedure is entirely analogous to that used in the analysis of the five-dimensional flows 
[TTJ |2H] , where the corresponding differential equations may be reduced to hypergeometric 
equations. In those cases, however, the regularity condition was imposed at a curvature 
singularity. 

A first inspection shows that (jB.l|) has two singular points (zero and infinity) with "s- 
rank {1; 3}" , in the language of [H]. This means that it may be obtained from a Fuchsian 
differential equation with four regular singularities by making three of these singularities 
coalesce at infinity. 15 This equation is known as the biconfluent Heun equation 021- With 

15 For comparison, the confluent hypergeometric equation has s-rank {1; 2}. 
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the change of variables 



x = e -^ u2 Y{u) , u = VPs , (B.3) 

equation (|B.1|) is mapped into a standard form 

uY" +(l + a- pu-2u 2 )Y' + {{j - 2 - a)u - \{5 + (3{l + a)))Y = 0, (B.4) 

with (a, P, 7, 5) = (a, 0, 0, 4P 3 / 2 ). The solution regular at s = is commonly denoted as 
N(a, (3, 7, 5; u). For /3 — 8 — 0, it reduces to a hypergeometric function; equation (jB.l|) is 
a different (but also very special) case. 

Despite considerable effort, see e.g. [32] and references therein, the Heun equation 
is still far less understood than the hypergeometric equation, which is the analogous 
fluctuation equation in the previously studied flows. For the purposes of computing 
correlation functions, there is even an additional technical complication here that comes 
from the fact that the relevant expansion (jB.2J) (i.e. the UV boundary of the flow) is 
around a generic regular point (s = 1), rather than around a singular point as in the 
higher- dimensional examples fTJ |2E]- Expanding around a singular point considerably 
simplifies the resulting expressions; the coefficient analogous to ^ in (jB.2J) for the case of 
the hypergeometric equation is usually denoted as ip and is simply expressed in terms of 
T-functions. In all, it is a tall order to solve (jB.l|) . but we find that we can extract the 
important information analytically, and use numerics to check. 

Following [12], we denote by H + (a, P,j,5;u) the unique solution of (jB.4|) that is 
regular as s — > oo. The relevant coefficient in (jB.2|) is then given by 

(B.5) 

Analytic expressions for H + (a, j3, 7, 5; u) may be constructed along the lines of jlH] as an 
infinite chain of sums over Pochhammer symbols (x) n = T(x + n)/T(x). 16 Rather than 
constructing these series, we will derive analytical results only for the IR asymptotics in 
P, and use the numerical solution of ()B.5jl for other purposes. 

To this end, we consider the following change of variables 

X = v- a/i Z(v), v = \P s 2 , (B.6) 

16 Unfortunately, the earlier results of 01] which seemingly give N(a, 0, 0, S) as a series in hypergeometric 
functions are incorrect; the double sum in eq. (2.8) in that paper does not, in fact, factor into (2.10); a 
relation which is recursively used in the construction. 



MP) 



-p 



o 



log H + (a,0,0,4P 3/2 ; VPs) 
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which transforms (jB.l|) into 



A a/4 Z = v 2 Z" + vZ'-(v 2 + ^Z = P 3/2 ^Z. (B.7) 

where prime now denotes a derivative with respect to v. This equation may thus be 
resolved into a series of inhomogeneous Bessel equations 

oo 

Z = p3r/2 2 ~ r/2 Z r > A Q/4 Z r = v 1/2 Z r ^. (B.8) 

r=0 

which can successively be integrated in terms of Lommel functions, demanding regularity 
in the interior v = oo. Note that this expansion is, in particular, compatible with the IR 
(small P) asymptotics. For instance, for a = —2, we explicitly find the first terms in this 
series as 



Z = y/2/nKi(v) = v~2e~ v , Z 1 = \^rv^e v evic(V2nj) , (B.9) 

with the complementary error function erfc(f). After some computation, this gives rise 
to the small P asymptotics 



= - P - P 2 + 2^P 5/2 -AP 3 + ... . (B.10) 

v=P/2 



For a = 0, one finds similarly the lowest term Zq = Kq(v), and thus the small p asymp- 
totics 

*° - C + lo 2 g(P/4) + — (BU) 

with Euler's constant C. As we have shown in the main text, the two-point correlation 
functions of the inert scalars in the 4 are proportional to the inverse of \I/o; the log-term 
then describes the standard behavior of a dimension A = 1 operator in the IR, cf. table UJ 
An analysis of the large P asymptotics of the gives 

m_ 2 = _v/3P + ! + ..., ^ = _V3P- i + (B.12) 

To obtain a numerical expression for the important ratio we can use the following 
simple prescription. Consider the differential equation (jB.4J) and numerically compute 
two of its solutions with two different sets of initial conditions given at u = yP (the 
regular point s = 1) as 

n(Vp) = o, y;(Vp) = -Vp, 

Y 2 (VP) = 1 , Y^(VP) = . (B.13) 
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Since there exists a unique solution regular as u — > oo (called H + {u) above), the ratio 
Y2/Y1 will tend to a constant in this limit, and this constant may be determined numeri- 
cally. From (jB.5|) . we then have the relevant coefficient ty a as 

* a = -P + lim J . (B.14) 

The result for a = —2, is plotted in figures 3 and 4 together with the first few terms 
of the exact asymptotics, given in (jB.10j) - (jB.12|) . Agreement in these regions is already 
quite good for including just a few terms. 




-0 . e 
-0 . 8 



-1 
-1.2 

Figure 3: Small P (IR) asymptotics 
of (straight), and \I>o (dashed). 
The dotted lines correspond to the 
first terms of the exact asymptotics 

(EE!), cnn. 




1 p 



Figure 4: Large P (UV) asymp- 
totics of l/^_2 (straight), and l/* 
(dashed). The dotted lines corre- 
spond to the first terms of the exact 
asymptotics (|B.12|) . The horizontal 
axis is 1/-P. 
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